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Summary. A simple exactly solvable kinetic model for the non-linear inelastic hard 
sphere Boltzmann equation is used to explore the relevance of hydrodynamics for a 
granular gas. The equation predicts a non-trivial homogeneous cooling state (HCS), 
including algebraic decay at large velocities. The linearized kinetic equation for small 
perturbations about the HCS is solved exactly. It is shown that the hydrodynamic 
excitations exist in this linearized dynamics, and their detailed form is shown to 
agree with results from the Chapman-Enskog method up through Navier-Stokes 
order. The existence of the hydrodynamic modes at short wavelengths, far beyond 
Navier-Stokes order is demonstrated as well. Finally, the precise sense in which the 
hydrodynamic excitations dominate the dynamics at long times is described. 1 



1 Introduction 

The applicability of a hydrodynamic description for a granular gas has been 
an issue of considerable interest in recent years P 03 El ■ The concepts in- 
volved and the precise questions to be asked have been clarified to a large 
extent in recent works 0], and partial answers to these questions have 
been given. The setting for all of these discussions has been the inelastic hard 
sphere Boltzmann equation. However, direct demonstration of the transition 
to a hydrodynamic description is made difficult by the analytic complexity 
of the Boltzmann collision operator. The objective here is to provide such an 
illustration by addressing these issues in the more limited context of a simple 
and exactly solvable model kinetic equation for a granular gas (SJ This 
model has the primary characteristics of the Boltzmann kinetic equation: dis- 
sipation of energy, a non-trivial homogeneous cooling state (HCS), and exact 

1 This research was presented at the "Modelling and Numerics of Kinetic Dissipa- 
tive Systems" workshop, Lipari, Italy, June 2004. An abbreviated version of the 
present work has been submitted for publication in the conference proceedings 
(L. Pareschi, G. Russo, G. Toscani Eds., Nova Science, New York.) 
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macroscopic balance equations for the hydrodynamic fields. Hence it is a good 
starting point for a demonstration of principle in the context of granular gases. 

A closed set of equations for the density, flow velocity and the granular 
temperature can be obtained from both the Boltzmann kinetic theory and 
the model kinetic theory considered here by the Chapman-Enskog procedure 
The drawback of this procedure is its practical limitation to states with 
small spatial gradients, and the assumption of a special form for the solution 
(a "normal" solution 0). There remain questions about the space and time 
scales on which such a solution can be expected. Also, the lack of energy 
conservation has led some to question whether the local energy or temperature 
should properly be considered a "slow mode" of the system that can dominate 
at long times. 

A direct probe of these questions is possible by characterizing all the ex- 
citations in the linearized dynamics of the system for small spatial perturba- 
tions of a homogeneous state, and looking for a time scale separation between 
the hydrodynamic excitations and the other microscopic excitations [2j. This 
prescription is difficult to carry out for the hard-sphere Boltzmann equation 
because of the complicated collision operator involved, even at the linearized 
level. However, this linear problem can be solved exactly for the model under 
consideration here. It is found that 1) hydrodynamic excitations exist and 
agree with those obtained from the Chapman-Enskog method under the same 
conditions of long wavelengths, 2) the hydrodynamic excitations exist over a 
range of shorter wavelengths far beyond the restrictions of the Navier-Stokes 
approximation, and 3) the other dynamics (called "microscopic" excitations) 
typically decay in a few collision times so that the hydrodynamic excitations 
dominate at long times. Some qualifications of 3) are noted. These conclusions 
hold even for conditions of strong dissipation. 

The paper is organized as follows. First, the definition of the kinetic model 
is recalled and its solution for spatially homogeneous states is given. It is 
noted that a special HCS solution results for a wide class of initial con- 
ditions after a few collisions times. Next, weakly inhomogeneous states for 
small perturbations about this HCS are considered. First, the hydrodynamic 
description obtained from the Chapman-Enskog results is recalled in order 
to define "hydrodynamic" excitations. Next, the kinetic equation is linearized 
and solved exactly by Fourier-Laplace transformation. The analytic properties 
in the complex plane determine all possible excitations, and it is shown that 
the hydrodynamic excitations are among these. Finally, the general solution 
is studied in the long wavelength limit where hydrodynamic and microscopic 
excitations can be displayed more explicitly. Conditions for the dominance of 
the hydrodynamic excitations at long times is considered in some detail. The 
results and remaining questions are summarized in the last section. 



Hydrodynamics for a granular gas 3 



2 Definition of the kinetic model 

A granular gas at low density is found to be well described by a system of 
N smooth hard spheres undergoing inelastic collisions parameterized by a 
coefficient of restitution a. The number density for particles with position r 
and velocity v at time t is determined from the Boltzmann equation [SJ |S] 

|+v-v)/(r,v,t) = J(r,v|/(t)). (1) 

The detailed form of the collision operator J (r, v|/(t)) is not required for the 
purposes here, only some of its properties. Among these are its moments with 
respect to 1, v, and v 2 




■7(r,v|/(t))= | |. (2) 



-§nTC, 



Here n, u, and T are the local density, flow velocity, and temperature at time 
t, obtained from the distribution function according to 

n(r, t) 

n(r,t)u(r,f) I = / dv I v I /(r,v,i). (3) 




]n(r,t)T(r,t) / 

Finally, £ is the "cooling rate" (see Ea. i|ll|) below to justify the terminology) 



C(r,t) = (l-a 2 ) 2An{r t)T{ ^ t) ] J dv ig 3 /(r,v,t,)/(r )Vl ,t). (4) 

The two zeros in Eq.@ correspond to the conservation of mass and momen- 
tum respectively. The source in the third term implies that the energy is not 
conserved and the system is "cooling" (note that it vanishes in the elastic 
limit of a — > 1). The functions n, u , and T are the hydrodynamic fields, 
and the properties J5J) together with the Boltzmann equation give the exact 
macroscopic balance equations relating these fields. These balance equations 
are the essential starting point for the formulation of a hydrodynamic descrip- 
tion for the system. Consequently, the properties should be preserved by 
any model kinetic theory intended to explore hydrodynamics. 

The kinetic model considered here is a simple extension of the BGK model 
for gases with elastic collisions 0, to capture the properties of a granular gas. 
It is obtained by replacing the bilinear collision operator in the Boltzmann 
equation by |S] 

J (r,v|/(t)) --*(/-<?), (5) 

where v is a velocity independent parameter of the model and g is chosen to 
be a functional of / such that the moment conditions are satisfied. As in 
the BGK model g is taken to be a Gaussian 
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g(v,V,t\f)=A(v,t)e- B ^ v2 , (6) 

where V = v — u(r,t) is the peculiar velocity. A and B are chosen so as to 
enforce the moment conditions above. They are found to be 

A(r,t)=n(r,t)(^) 3/2 , B(r,t) = m (7) 

where m is the mass. In the elastic limit g becomes the local Maxwellian; 
otherwise, the effective temperature is changed to account for the cooling 
implied by Eq. J2J). 

It remains to choose the functions v(r, t) and ((r,t). In principle, these 
are specific functionals of / in the Boltzmann equation. Here, they are taken 
to depend on / only through the temperature and density. The cooling rate is 
chosen to be the same as that obtained from the Boltzmann equation by using 
a local Maxwellian for / in Eq.Q) The collision frequency v is chosen so 
as to fit one of the Chapman-Enskog transport coefficients of the model to 
that of the corresponding result from the Boltzmann equation, in this case 
the shear viscosity rj (for details see [§])■ These choices give 

C(r,i) = ^(l-a 2 )^ (r,t), u(i, t) = - 1(1 - a) 2 ^o(r, i), (8) 
where ^(r, t) is an average local collision frequency 

v (r,t) = — n (r,t W V ; . 9 
5 V m 

This completes the definition of the model kinetic equation 

f^+v-v) / = -*/(/-<?). (10) 

This is deceptively simple, since the right side is a highly nonlinear functional 
of / through its dependence on the hydrodynamics fields in v and g.. 



3 Spatially homogeneous states 

As in the case of the Boltzmann equation, there exists no spatially homoge- 
neous steady solution to 110|) for the isolated system. This can be seen by 
taking the second velocity moment of the kinetic equation to get 

r- 1 a t T = -c. (ii) 

This implies that one of the moments of the distribution function is always 
time dependent and hence there is no steady solution. Instead, there is a 
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special scaling solution, the HCS, which is approached in a few collision times 
by most spatially homogeneous initial conditions. 

Consider a general homogeneous initial distribution and look for solutions 
to the model kinetic equation in the dimensionless form 

/ (v, t) = nv^J* (v* ,«), v* = v/v hcs (t) (12) 

Here the velocity scaling is generated by the time dependent thermal velocity 
and s is an average collision number 



v hcs (t) = y/2T(t)/m, ds = is (t)dt. (13) 
The dimensionless form for the model kinetic equation becomes 

d s + ^C*V V , • v* + v*^j f* (v*, a) = v*g* (v* | f ) . (14) 

This can be integrated directly for initial homogeneous states to get the gen- 
eral solution (see Appendix A and |§]) 



/*(v*,s) = e-(^" +l/ ') s f*(e-^* s v*,0) + is* (tt (\ - i 



3/2 



ds'e-&* +,/ *>' exp ^ ~ j i (15) 

with the dimensionless constants ^* = i^/^o, and = C/^o- Since the state 
is homogeneous the flow velocity has been taken to vanish, u = 0, by making 
the appropriate Galilean transformation. 

The constant v* is of order unity, so the domain of integration in l|15l) is 
exponentially bounded for s > 1 and for large s the integral becomes inde- 
pendent of s. The first term vanishes exponentially fast for initial conditions 
uniformly bounded in the sense /*(e _ 2? s v*,0) < C. Therefore an s— inde- 
pendent HCS solution is obtained 

"' (16) 

In terms of the dimensionless variables chosen, the HCS behaves as a sta- 
tionary state. It is universal in the same sense as the Maxwellian for elastic 
collisions, since most initial homogeneous states evolve to it after a few colli- 
sions. In the following it will be chosen as the reference state to study small 
spatial perturbations of a homogeneous state. Clearly, the results obtained 
will be similar for other homogeneous reference states, but with an additional 
transient period as the reference state itself approaches the HCS. 
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4 Hydrodynamic excitations 

Hydrodynamic excitations are most easily denned in the long wavelength limit 
where the Navier-Stokes hydrodynamic equations apply. Consider inhomoge- 
neous states characterized by smooth spatial and temporal variations in the 
density, temperature and flow velocity, and for which the deviations of the 
hydrodynamic fields from their values in the HCS are small (linear hydrody- 
namics). The detailed form of these equations and the transport coefficients 
occurring in them have been obtained from the kinetic equation (Boltzmann 
and kinetic model) by the Chapman-Enskog method ^] and are found to 
be 

d t Sn + n hcs V ■ 5u = 0, (17) 

d t Sui H — (n hcs diST + T hcs di5n) 

n hcs m 

' n (djSui + diSuj - h tJ W ■ Su ) ) = 0, (18) 



nh cs m \ \ 3 



d t 5T + \T hcs SJ ■ Su + ICkcsST + Ih^Sn 



2 'nV 2 5T + nV 2 6n) = 0, (19) 



3n h 

where 8x = x — x^ca denotes the deviation of a hydrodynamic field from its 
value in the HCS. The transport coefficients 77, K,and fi occurring in the above 
equations have been evaluated in reference jS] . 

The linear Navier-Stokes equations are complicated only by time depen- 
dent coefficients due to the cooling reference state. This time dependence 
can be eliminated by introducing the dimensionless variables, Sn* — 5n/rihcsi 
5T* = ST/Thcs: and Su* = Su/v^cs- In addition the time is transformed 
to the collision number s in 1)13(1 and the dimensionless space variable is 
r* = vvq/vhcs- Finally, since the equations are linear, it is sufficient to consider 
a single Fourier component 

y:(r, S ) = e ikr *y:(k, S ) 

where the hydrodynamic fields are chosen to be 

S£(k,s) «-> (*n*(k,a),<yT*(k,a),<5u*(k,s).k,<yu*(k,a)-ei,<yu*(k,«)-e 2 ) . 

(20) 

The components of the flow field are taken to be those along the orthonor- 
mal triplet {k,ei,e 2 }. The dimensionless linear Navier-Stokes equations then 
become 

S s y*(k,s)+/C a ^(k, S ) = 0, (21) 
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with 



K 



( ik \ 

C* + f M*fc 2 £■ + |K*fc 2 fife o o 

if " j| -^ + §77*fc 2 

-^ + i^*fc 2 



2' 

V -^- + ±T]*k 2 J 

(22) 

The transport coefficients are scaled with respect to their values in the elastic 
limit, k* = k/ko, ?y* = r)/r]o, H* = (MhcsThcs/ '^0; where the elastic limit values 
are n = 157y /4m and 770 = 5(mT^ cs ) 1 / 2 /16o' 2 7r 1 / 2 . The dependence of 77, K, 
and n on the restitution coefficient a has been evaluated in reference [S]. 
The initial value problem can be solved directly with the solution 

y* a (Ks) = J2 C ^ (k)S - (23) 



Here C a p are constants and zp(k) are the eigenvalues of K. (k) determined 
from 

det(z(fc)l+/C(fe)) =0. (24) 

As the Navier-Stokes equations are valid only up through order k 2 the solu- 
tions to this equation are relevant only to the same order 

1,2 C* , ( 2 5 A 2 C* / 1 , 2 A - 



^ ^ h ' 2 + ^ - i K V * ' 2 ~ lac + * 

For the discussion below it is noted that in particular, at k = 0, these become 

The hydrodynamic excitations can now be defined precisely. They are the 
excitations of the form (|23|) with zp (k) defined to be the frequencies that 
are continuously connected to z a (0) defined in l|26|l as k goes to zero. When 
such modes can be identified it is said that hydrodynamic excitations exist. 
In that case, an internal consistency check is the agreement of those modes 
with the Chapman-Enskog results H25|) at order k 2 . The modes 24 (fc) = z$ (k) 
are associated with the components of the flow field orthogonal to the k. 
These are called transverse shear modes. The other three modes are called 
longitudinal modes. Since K. is block diagonal the transverse modes decouple 
from the other modes and hence can be treated more easily and directly. 
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5 Linear kinetic theory 

In the rest of the paper, the dimensionless variables introduced in the previ- 
ous sections are used throughout and the asterisk is deleted for simplicity of 
notation. To describe small perturbations of the distribution function about 
the HCS solution, define Sf by 

/ = hcs + Sf. (27) 

Substitution of this form into the kinetic equation l|10|) and retaining only 
linear order terms in Sf defines the linear kinetic theory. The model collision 
operator of H10|) depends on the distribution function explicitly and implicitly 
through the hydrodynamic fields in v and g. The linearized kinetic equation 
is identified in Appendix B and is found to be 

d a + vV T + v+ ^CV V • v^jSf = M a y a (28) 

where 

/ (CV V • v + v) f hcs \ 

M Q (v) = ±V V • v (|C - (K v v ' v + u)) f hcs . (29) 

V (-5C-(|CVv'V + I/))Vv/hc. / 

The hydrodynamic fields y a are defined in l|20(l . The kinetic equation can be 
solved by a Fourier-Laplace transformation 

6f(k, v, z) = J dse- zs J dre lkr 6f(r, v, s). (30) 

The formal solution is 

S~f(k, v, z) = KSf(k, v, 0) + g a (k, z)KM a (v) (31) 
where the resolvent operator 1Z is 

K(k,z) = ^z + ^ + ik-v+icVv-vj (32) 

This is an implicit solution for <5/(k, v, z) because the hydrodynamic fields 
{y a } are linear functional of Sf 

y a (. k ,z) = J dv<j) a (v)<5/(k,v,z). (33) 

where 

1, ( ^v 2 -l) ,vk,v-ei,ve 2 J (34) 



3 

The {y a } can be determined self-consistently using (|3*T)l in (|3^|) to get 
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$ a (k,z) = (l-Q(k,z))- 1 pIp(k,z) (35) 

where 

I a = / dv^(v)^/(k,v,0), Q a/3 = J dv<t> a {v)1lMp. (36) 
Finally, substitution of ((35(1 in 131(1 gives the desired solution 

*/(k,v,z)=7W/(k,v,0)+7JM o (v)(l-Q)^I / s(k,«). (37) 

It remains to make explicit the action of the resolvent operator in ((3711 . 
This is determined in Appendix A with the result, for an arbitrary function 
of the velocity X (v) , 

KX (v) = J dse~( z+u+ i c ) s exp Q (V^ Cs - l) ik • v) X (e^v) , 

(38) 

In this way the general linear solution (|37() is reduced to quadratures. / and 
Q in Eg. 136(1 are evaluated in Appendix C. 

5.1 Transverse modes 

The analysis for transverse shear mode excitations is simplest and illustrates 
the more general case. Consider an initial perturbation in the transverse flow 
velocity only 

5f(k, v, 0) = -e x • V v f hcs (v) Su t (k, 0) u t = u-ej (39) 

This perturbation results from a small change in the velocity of the HCS along 
the transverse direction e\. It couples only to uei (due to symmetry) and 
Eq. 1(35(1 for the response of that held reduces to 

TYk z) 

8u t {k,z) = K ' I gfi t (k,0). (40) 

with /(k, z) = -?4(k, z)/Sut(k, 0). The functions J(k, z)and Q^Ck, z) are eval- 
uated using the property ( I38|) with the results 

/*°° r°° 
I{k, z) = J dse-( z+iy -^) s v J ds'e~ vs ' exp [-e^' a{s)k 2 J , (41) 

where 

Similarly, Qa{k, z)is found to be 
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Qu(k,z)=[v-tk^)I(k,z) (43) 



Consider first the long wavelength limit k — > 



z — 



1(0, z) = l_Q(0,z) = — * (44) 

z + ^ — | z + — | 

«fit(k, z) -» — ^-afitfr 0). (45) 

Comparison with shows that this is the expected hydrodynamic pole, 
confirming the existence of hydrodynamics in the kinetic theory for response 
to a transverse perturbation. 

The hydrodynamic pole in the long wavelength limit arises from the zero 
of 1 — ^44(0, z(0)). Its continuation to z(k) at finite k is therefore defined to 
be the solution to 

1-Q 4 4(k,z(fc))=0, (46) 

with the condition z(0) = C/2 (the poles of the function /(k, z) determine the 
microscopic excitations in <5u t (k, z) for finite k). It is easily shown that z(k) 
is a real function of k 2 . The integrals above are defined only for z > | — v. 
Therefore, the hydrodynamic pole can exist only for values of k such that 
z{k)>-(v-(). 

Figure 1 illustrates the numerical solution to (|46|l for a = 0.8, showing the 
hydrodynamic pole exists for all k < k m , where z (k m ) = | — v. For small k 
direct expansion of l|46() to order k 2 gives the solution 

*(*) - l<(<*) ~ \ri{a)k 2 . (47) 

This is the same form as that for the transverse hydrodynamic mode obtained 
from the Navier-Stokes equation in (|25|l . The shear viscosity obtained here is 

= ( ( \ l -T7( \v (48) 
(v(a) - ^Cw) 

which is the same as that found by the Chapman-Enskog method in |S] and 
hence is an independent verification of that method. Interestingly, Figure 1 
suggests that the Navier-Stokes approximation begins to fail for relatively 
small wavevectors. 



5.2 Longitudinal modes 

Return now to the general solution and consider initial conditions that 
generate the perturbations 6h(k., 0), <5T(k, 0), and <5u(k, 0) • k. The Laplace 
transform of these fields at later times is given by Ij35(l 
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Fig. 1. Transverse hydrodynamic mode 



~y a (k, z) = (l- Q(k, z))-^I (k, z), a, (3 = 1, 2, 3. 



(49) 



The integrals z)and Q a p(k, z) can be calculated exactly just as for the 
transverse perturbation above. As in that case, the poles of (1— Q(k, z))~ 1 gives 
the hydrodynamic excitations. In place of (|4(j[l . the longitudinal hydrodynamic 
modes are now defined by the three solutions to 



det(l - Q(k,z(k))) = 0. 
At k = 0, the matrix 1 — Q(k, z) is found to be 



(50) 



l-Q{Q,z) = 



\ 



2 + 1/ 

z+u-C ^ 







z-\-u z-\-v—Q 





(51) 



z-\-v- 



7 



The solutions to (|5U|I give the longitudinal modes in the long wavelength limit 



z(0) 



' 2' 2 



(52) 



Comparison with (|2(j|) shows that these are indeed the hydrodynamic modes. 

It is straight forward to extend the calculation of these modes by expanding 
Q(k, z) to order k 2 . The solutions to (|5t)|l are found to be the same as in 125|l 
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1,2 C , (2 5 \ C fl 2 > ,_ 2 



*(k) - (0 + ^, -5 + ^ - r ) f - ^ + ~ 3 V) k ' j . (53) 

The shear viscosity rj is again given by 1)48 [I. while the thermal conductivity 
k is 

2((Ka)-C(c*)) 2 +C 2 W) 
" (a) = 3(K«) " 2C(a))Ma) " ( > 

This agrees in detail with the expression obtained by the Chapman-Enskog 
procedure 

The results of this section complete the first objective of this work: demon- 
stration that the hydrodynamic excitations exist in the solution to the linear 
kinetic equation for small perturbations of the HCS. The form of the hy- 
drodynamic modes to Navier-Stokes order and dependence of the transport 
coefficients on the restitution coefficient is confirmed in detail, providing in- 
dependent support for the assumptions of the Chapman-Enskog method. In 
addition, Figure 1 shows that the hydrodynamic modes exist far beyond the 
validity domain for the Navier-Stokes approximation. For the example con- 
sidered here that domain appears to be quite restricted. 



6 Approach to hydrodynamics 

The second objective of this work is to verify that the hydrodynamic excita- 
tions are the slowest modes at long wavelengths and therefore dominate at 
long times. First, the general solution to the linearized kinetic equation 1371) 
is reconsidered to display all of the possible excitations. 

<5/(k, v, z) = TZSf(h, v, 0) + ftM a (v)(l - Q(k, z))-%Ip{k, z). (55) 

The dynamics arises from two types of contributions. The first is the non- 
hydrodynamic behavior generated by the resolvent operator in TZ5f, lZM a , 
and i/3(k, z) (see Eq. H38|| s |. The hydrodynamic exciations arise from the simple 
poles of (1 — Q(k, z))~p. Also, there is a convolution of the hydrodynamic and 
non-hydrodynamic behavior. All of these effects are present even for k — > 
0, where these expressions can be made more explicit. Hence the following 
analysis is first carried out at k = and then it is argued that the conclusions 
extend to k ^ by continuity. 

Using (|38|l and the results of the previous section at k = the general 
solution can be written exactly 

<5/(0,v, s ) =^ e ^(°) s XQ (v)^(0,0) + e -(^) s (*/(0,e-*<«v,0) 

a 

l^e z »(»> Xo ( e 4(^( ) , (56) 

a ) 
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The new hydrodynamic fields {y' a }in (|56|l are linear combinations of {y Q }and 
the functions Xa (v)are the corresponding linear combinations of M a (v) 



lies 



}• 



y' a = j<fo, 5n + <5u j , Xa = {{fhcs + V v • vfhcs) , -V v ■ vf hcs , V v // 

(57) 

The first term on the right side of l|56[) is the contribution from hydrodynamic 
excitations. It could have been constructed directly from the solution to the 
linearized kinetic equation as an eigenvalue problem. Then Xa( v ) would be 
the hydrodynamic eigenfunctions. This has been done in [3] and the above Xa 
are indeed the eigenfunctions found there. 

All other excitations in the system are contained in the second term of (|56|l . 
Nominally, this appears to decay as exp(— (§C + i/ ) s )- This is misleading, 
however, due to the additional time dependence in the functions of velocity. 
Consider first the moments giving the hydrodynamic fields y' a (k, s) 

y' a (k,s) = Jdvxl> a (v)6fQi,v, S ), i> a = jl, X -v 2 + ~ v| . (58) 

Then l|55) l gives 

^(0, s) = e^°>y'JO, 0) + e-"' J d^ a (e^'v) {sf(0, v, 0) 

+ Y,e Z ' mS Xp(v)y' P (O,0)) (59) 

a J 

Since the functions ip a (v) are polynomials of maximum degree 2, the second 
term decays as exp(— [y — Q) s). The fastest decaying hydrodynamic mode is 
that with 2:2(0) = — C/2- Thus the hydrodynamic modes are bounded away 
from all of the microscopic modes and the fields y' a (0, s) obey the hydrody- 
namic equations for long times, if v > 3£/2. This must be the case since this 
condition characterizes energy loss, according to the moment conditions on 
the collision integral For the model considered, 3£/2 =5 (l — a 2 ) /8 and 
v = 1 — (1 — a) 2 /4, so the condition holds even at a = 0. 

Now consider other properties obtained as averages over the distribution 
function. In general, for a property determined from A(v) 

(A;s) = (A(v)) hcs +J2e z « i0)s y' a (0) [ dvA(v) Xa (v) 

+e- us J dvA(e^ s v) ^/(0, v, 0) + £ e z ° (o)s Xa(v)^(0)^ (60) 

For any bounded function |A(v)| < G the last term decays as e~ vs , and 
since v > Q/2 the hydrodynamic contribution dominates for long enough 
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times. However, a peculiarity occurs for unbounded functions of sufficiently 
high degree. Suppose A(v) ~ v p for large v. Then the last integral will have 
a contribution behaving as exp(— (y — |£) s). The maximum value of p is 
restricted by the condition that the average in the reference state, (^(v))^, 
should be finite. This condition gives p < Iv/C,, so the exponent for the second 
term remains negative. Nevertheless, the decay rate can be made as small 
as desired within this restriction. In particular, it can be chosen to give a 
decay that is slower than the fastest hydrodynamic modes at £/2. For such 
properties, the hydrodynamic excitations contribute but they do not dominate 
at long times. 

The conclusion of this last paragraph suggests new limitations on hydro- 
dynamics for granular gases. There are two ways to avoid such a conclusion. 
The first is to assume that such properties that behave as A(v) <~ v p with 
{y — |£) < C/2 are unphysical or uninteresting. For example, at a = 0.9 this 
gives p > 11. A more precise avoidance of the problem is to formulate the 
linear kinetic theory problem in an appropriate Hilbcrt space whose scalar 
product is defined by 

(a, b) = J dvf hcs (v) a* (v)6(v). 

The solutions of interest are then Sf = fh cs A., where A must lie in the Hilbcrt 
space. Similarly, the properties whose averages are considered must also be 
in that Hilbert space. Then square integrability changes the above restriction 
on p to 2p < 2v/(. This restores the dominance of hydrodynamics at long 
times for all admissible physical properties. The choice of a Hilbert space may 
seem artificial since normalization of the distribution function requires only 
existence in L\. However, for normal gases the Hilbert space has been justified 
on the basis of the existence of the H functional, or entropy, when linearized 
about the equilibrium state. Perhaps there is a corresponding functional for 
granular gases leading to the same conclusion. 

7 Summary and conclusions 

Granular gases are significantly different from normal gases. Energy is not 
conserved. So, the homogeneous states for an isolated system are time depen- 
dent and the velocity distribution is not Maxwellian. Such states are unstable 
with respect to small spatial perturbations. These peculiarities have often 
been cited as reasons to question the relevance of hydrodynamics for granular 
gases, or at least to expect severe limitations on the applicability of such a 
description. The model kinetic equation considered here supports all of these 
peculiarities and therefore provides a good testing ground for such concerns. 
Previously, the model has been used to derive hydrodynamic equations by the 
Chapman-Enskog method for states with small spatial gradients. The analysis 
here is complementary, in the sense that it leads to the same results under 
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those restricted conditions but describes as well the case of large gradients 
and the time scale leading up to the hydrodynamic stage. 

The problem posed here is to solve exactly the linear kinetic equation for 
small spatial perturbations of the homogeneous reference state (HCS). This 
has been accomplished by first solving the equation for the HCS and then con- 
structing a solution to the linear equation valid for arbitrary space and time 
scales. It has been shown that hydrodynamic excitations exist in this gen- 
eral solution by identifying them explicitly at long wavelengths and assuming 
an analytic connection to shorter wavelengths. In this way the Navier-Stokes 
hydrodynamic excitations are demonstrated explicitly. The existence of hy- 
drodynamic excitations at shorter wavelengths was illustrated for the special 
case of the shear mode. In the last section above, it has been shown that 
these excitations dominate the other dynamics for the hydrodynamic fields 
and other average properties after a few collisions. In summary, it appears 
from the results of this model kinetic equation that the existence and condi- 
tions for a hydrodynamic description for granular gases is similar to that for 
normal gases. This conclusion is restricted to the case of small perturbations 
of an isolated gas, and further study is required for the more general and more 
interesting cases of boundary driven systems and external forces. 
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A Solution of the kinetic equation 

The formal solution to the model kinetic equation 1|14|) is 



The action of the exponentials in (|61|l can be determined as follows. Let us be- 
gin by considering the more general case that includes spatial inhomogeneities. 
Define a function X(v* , fc, s)by 



/*(«*, s) = e -(i C{3+v '- v -' )+iy *> f*(v* ,0) 




(61) 



X(v*,k, s) = e -(K*(3+v*.V v .)+^+ik.v) SX ^ 



(62) 



which then obeys the equation 




(63) 
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Next introduce 

X(v*, k, s) = e-^* sv *- v -"X(v*,k, s) (64) 
so that X(v* , s)obeys the equation 

d s + |c* + e^* sv *- v -* [y* + ik • v*) e -K*«v*-V v , X = 0. (65) 

From the identity 

e iC*™*-v v , F ( v *} e -±Csv*-v v , = F(e* c * a v*} (66) 

this equation becomes 

9 s + ^C* + ^*+e^ C * s ik-v*^X = (67) 

This can be integrated directly and inserted in l|64|l to give 

X(v*,s) = e-i c * s e-^* sv * v -*e-^ s exp (- J ds'e^" s 'ik • v*^) X (v*) 

= e-i Cs e-'* s exp(- ds' e -^ c s ' ik ■ X (e~^ c s v*^ (68) 

In particular, if k — 0,thc formal solution to the kinetic equation for a homo- 
geneous system becomes 

/*(v*,s) = e~ ! i Cs e- / * s f*(e-^' s v*,0) 

+ f ds'e-%< : * s 'e- / * s 'p*g*(e-^* s 'v*,s- a'), (69) 
Jo 

B Linearization of the model kinetic operator 

The collision operator Ea. (|10f> to linear order in the small deviation 6f be- 
comes 

J[f] -> -Vhcs (fhcs - 9hcs) - Vhcs (Sf - Sg) - 8v (f hcs - g hcs ) (70) 
Now, since v and g are normal 

x x dg dg hcs 

Og = OVa-K \hcs= d Va-^ |v=v, 

5v = Sy a -^- \ hcs = Sy a ^ hcs | v =v (71) 
oy a oy a 

where y^s are the hydrodynamic fields as before. Note that the derivatives 
must be taken using the HCS as a function of V = v — u, and then setting 
u = 0. The kinetic equation becomes 
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dy a 



d 1 df 

(d t + v • V+^cs) Sf = -Sy a \ -g— [v (f - g)} j | hcs +v hcs 8y a 

I 

dfhc 



h cs 



Sy a 

+Vhc S 5y a 



d 1 

q oChcsVv • (Vfhcs) 

dy a 2 



V=v 



(72) 



Use has been made of the fact that the HCS obeys the equation 

TrOtcs-Vv ' (V// lcs ) = V (fhcs — 9hcs) ■ 



(73) 



in the second equality. More explicitly, carrying out the derivatives using the 
variables n, T,and Ui 



d 

(d t + v • V+Vhcs) Sf = Sn- 
ort 



^OicsVv • (Vfhcs) 



V=v 



-Su • V u 



ic (0) V v • (Vfhcs) 

k (0) V V ■ (Vfhcs) 



v= 



V=v 



i j dfhcs j- dfhcs , crr^fhc: 

+Vhcs Sn— h Siii— 1- ST— — 

on ouj dl 



(74) 



Using the normal form of the local HCS these derivatives can be evaluated to 
give 



(d t + V • V+Vhcs) Sf = (Ch cs V v • vf hcs ) 



Sri 



rihcs 

i St 

"tOics (V v ' Vfhcs - V v • VV V • v/ft c5 ) — 

4 -I has 

i Sti 

'T^Chcs (y vfhcs + V v ' vV v // lcs ) • <5u + Vhcsfhcs 



V hcs~ f ST 

~ Vv • VJhcsTf VhcsVvfhcs ■ OU 

^ 1 hcs 



(75) 



This equation in scaled variables becomes 



d s + ^C*V V . • v* + v* ■ V r . + v* j -)./' 1 : : 1/ , ( v ' i ,/, , i r 1 . .S) I 71, i 



where M*(v*) and y*(r*,s) are as defined in Eq.J2H) and Eq. (120)1, respec- 
tively. This gives the linearized dynamics of the system for small deviations 
from the HCS. 



18 Aparna Baskaran and James W. Dufty 

C Evaluation of / and Q 

Here, the action of the resolvent operator defined in Eq. (|32[) of the text is 
made explicit. In particular, the functions 

I a (k,z) = f dvct> a (v)1Z5f(k,v,0), Q Q/3 (k,z) = J dv<f> a (v)KMp. (77) 

are evaluated. Note that scaled variables are used and the superscript asterisks 
have been suppressed. First, define a Fourier transform with respect to the 
velocity v of the form 

Y(X) = { dve lX - v Y(v). (78) 



Then, Eg. (|77|) above can be written as 

I a (k,z) = (27r)~ 3 / d\s} t (k,X,0) / dv</) Q ^e~ jA v , (79) 



and 

Q a p{k,z) = (2tt) _3 J dXMp (A) J dv^Tle^^ (80) 

Hence in general it is sufficient to describe the action of the resolvant on the 
Fourier mode function e~ zX v . 

It follows from (|62() and l|68[> Appendix A above that 

poo 

Jo 

= / dse- 2fl e-^ fl e-" s exp(-ik(s) -v^-* 6 " 5 A ' v , (81) 
Jo 

k(s) = f ds'e-i^'k. (82) 
Jo 



with 



It follows that 



J d^ a Tle- lX - v = J™ dse- (z+v)s e-^ s J dv^ a (v) e -<A-e i«'v e -<k(.}.v 

= J™ dse~ {z+v)s J dv0 Q (e^ s v) e *(H-s)-V-v (83) 
Putting in the forms of the {<fi a } given in Eci . (|34l ) leads to 

/>00 

dv^Ke- ,A v = / dse^ z+v)s (2ir) 3 6 {k{- s) - A) (84) 



dvfolZe 



-i\-v 



(°° dse~ {z+v)s - l) i^f 5 (k(-s) - A) 

(85) 
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dv(f> 3 Ke- lX - v = ik ■ V A / dse-( z+ "-3C) S ( 2?r ) 3 5 ( k (_ s ) _ A ) (86) 



V A / dse-(^-iO s (27r) 3 ,5(k(- S )-A) (87) 



This allows evaluation of J Q (k, z) 



/i(k,*)= / d S e-( 2+t, *^/ f (k,k(- S ),0), 



7 2 (k,z) = jf d S e-^*) s (-- e <*s V £ ( _ s) - lj <5/ t (k,k(- fl ),0), 

/"OG A 

/ 3 (k,z) = / d S e-(^-^*) s ik.V k( _ s) 5/ t (k,k(- S ),0), 
Jo 

h, 5 (k,z) = / d Se -( z+ "*-5C*)^ 12 . Vk( ^ )( 5/ t (k,k(- s ),0), 
Jo 

The elements Q^^k, z) are evaluated in a similar way 

/>oo 

Qip(k,z)= dse- {z+v)s Mi3 (k(-s)) 



Q 2/3 (k, z) = d se -^> s (-HgC-v^ - l) M (k(- S )) 

poo 

g 3/3 (k, *) - / dse-( z+v -i£) s ik ■ V k( _ s) M^ (k(-s)) 
Jo 

/■oo 

Q 4 ,5/3( k , ^) = / dse-( z+v -^) s ie h2 ■ V u( _ s) M^ (k(-s)) 
Jo 

The (A) are obtained directly from the definition of {Mp} in Ea. l|29(l 



(89) 
(90) 
(91) 

(92) 

(93) 
(94) 
(95) 



I (-CA ■ V A + v) f hcs (A) 

-±A • Va (|C ~ (~K A ' V * + ")) ^hcs (A) 

V (K (1 - |A ■ V A ) + 1/) iXhcs (A) 



(96) 



where 



Ac, (A) = / dve^f hcs (v) 



(97) 



dwe 



?'A-v 



"•»>(-) e V J 
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with the definition 



A(s') = — =A. (98) 



As an example, Q44(k, z) for the transverse shear mode is 
Q 44 (k, z) = J™ d S e-( z+v -H> (-ie, ■ V A Q( (l - \^ ■ V A ) - ^ (99) 

xiei • Xfhcs (A)) |A=k(- s ) 

(j \ rOO rOC 
v -tfJL\ I d8e-( z+ "-W a v ds'e- vs ' exp (-e^'a(s)fc 2 ) 

(100) 

where 



a(*)=^l-^] ( j . (101) 

This is the result given in the text. 
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